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Abstract 

In an interesting class of models, non-renormalizable terms of the superpotential are responsible 
for the spontanteous breaking of Peccei-Quinn (PQ) symmetry as well as the generation of the 
term. The flaton fields which break PQ symmetry are accompanied by fiatinos, and the lightest 
flatino (the LSP) can be the stable, while the decay of the lightest neutralino (the NLSP) might 
be visible at colliders with a low axion scale. We examine the cosmology of these models involving 
thermal inflation just after the PQ phase transition. The branching ratio of flatons into axions 
must be small so as not to interfere with nucleosynthesis, and flatons must not decay into the LSP 
or it will be over- abundant. We explore a simple model, with light flatons which can decay into Z 
or W bosons, or into a light Higgs {hP) plus a Z boson, to show that such features can be realized 
in a wide range of parameter space. The mass of the NLSP can be as low as (m/jO +m^)/2, with an 
axion scale of order 10^*^ GeV and a final reheat temperature typically of order 10 GeV. Then, the 
flatino LSP is a good dark matter candidate because the reheat temperature can be high enough 
to allow its production from the decay of the thermalized LSP, while low enough to prevent its 
overproduction from the decay of sfermions. 
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I. INTRODUCTION 



In supersymmetric extensions of the Standard Model, both the problem and the strong 
CP problem can be resolved through a simple extension of the Higgs sector, which im- 
plements spontaneously broken Peccei-Quinn (PQ) symmetry [1]. Suppose that the Higgs 
bilinear term in the minimal supersymmetric standard model (MSSM) superpotential takes 
the non-renormalizable form; 

W^.-term = h^H.H, (l) 

where Mp ~ 2.4 x 10^*^ GeV is the reduced Planck mass, p = J2ih ^ and PV^-term respects 
a PQ symmetry. Then, nonzero vacuum expectations values (vevs) can spontaneously 
break PQ symmetry, providing the axion solution to the strong CP problem [2], while 
generating the right order of magnitude for the /i parameter [3] ; 

f = A^. (2) 

In the case of a polynomial of the second order (p = 2), with /i ~ 1, this gives n ~ F^/Mp ~ 
100 GeV if the axion scale Fa ~ is of order 10^*^ GeV, its lowest feasible value. With a 
higher power p > 2, one would need a larger value of Fa to get the right size of 

From now on, we consider the case of supergravity with gravity-mediated supersymmetry 
breaking [4] . There are two ways of spontaneously breaking PQ symmetry. One is to intro- 
duce a renormalizable tree-level superpotential. In this case, the particles corresponding to 
the fields and their spin 1/2 partners, all have mass of order Fa except for a supermulti- 
plet comprising the axion (the pseudo-Goldstone boson of the PQ symmetry), the spin zero 
saxion, and the spin 1/2 axino. The mass of the axino in this case is generically of order 100 
GeV, but it can be as small as O(keV) [5, 6]. If it is the lightest supersymmetric particle 
(LSP), there follow some interesting consequences in cosmology and coUider physics. If the 
axino mass is ^ 10 GeV, it can be the cold dark matter of the universe if its (non-thermal) 
production comes mostly from the decay of the next-to-lightest supersymmetric particle 
(NLSP) [7]. For this to happen, of course, the primordial LSP relic density has to be diluted 
away, which would be the case when the reheat temperature after infiation is sufficiently 
low [8]. (In contrast, the axino mass is required to be less than a few keV if its primordial 
abundance is not washed out [9].) 

The alternative to a renormalizable superpotential is a non-renormalizable superpotential 
with soft supersymmetry breaking, leading to a fiaton model of of PQ symmetry breaking 
[10]-[22] . (A fiaton [16, 23] is a scalar field with vev much bigger than its mass, corresponding 
to a very fiat potential.) For simplicity, the superpotential is supposed to be dominated by 
a single term, 

^p^=jl^,m (3) 

where n — Y^iki > ^. If the soft mass-squared of one or more of the singlet fields (j)i is 
negative, e.g., through a radiative mechanism [15], there exists always a nontrivial global 
minimum. The (pi with nonzero vevs are then flatons. The axion scale is F.^ ^ ~ 
(moMp^^)^/""^, which is naturally of the right size. For instance, for n = 4, / ~ 1 and 
mo ~ 100 GeV, the vacuum expectation value is of order 10^° GeV, but it can be bigger if 
n is bigger. 
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In a flaton model, the masses of the particles corresponding to the fields (pi, and their 
spin 1/2 partners, all have mass of order 100 GeV, with the sole exception of the axion. 
The saxion and axino, defined as the superpartners of the axion, have no special significance 
and are in general not mass eigenstates. Instead, with N >2 fields breaking PQ symmetry, 
there are 2N — 1 spin zero fiaton particles, and N spin 1/2 fiatinos. These play the same 
role as, respectively, the saxion and the axino of renormalizable models. In particular, the 
lightest fiatino may be the LSP and hence a dark matter candidate. 

With either an axino or a fiatino LSP, the NLSP, which is typically the hghtest neu- 
tralino in the MSSM, decays to the LSP with the rate proportional to An interesting 

implication for collider experiments has been pointed out in Ref. [24]. Although the NLSP 
decays are very weak, and thus the corresponding decay lengths can be much larger than a 
collider size, there is an opportunity to observe these decays in future colliders with a large 
number of supersymmetric events so as to provide direct information about physics at very 
high energy scale. 

In this paper, we study the cosmology of models with a fiatino LSP, identifying various 
bounds on the parameters which are required for a viable cosmology. The paper may re- 
garded as a sequel to [22], in which all of the fiatinos were supposed to be unstable. We 
shall focus particularly on the case of a low axion scale. Fa ~ 10^° GeV, which increases the 
possibility of collider signatures. 

II. COSMOLOGY OF FLATINO LSP 

Flaton models lead to very interesting cosmological effects. On the reasonable assump- 
tion that the fiaton fields are nonzero in the early Universe, they lead to thermal infiation 
which can dilute away all the unwanted relics [16, 21]. In a fiaton model of PQ symme- 
try breaking, some restrictions have to be put on the model in order not to overproduce 
unwanted relics again from fiaton decay. As a (scalar) fiaton field (p has always the decay 
channel ip —>■ aa with a rate T^^aa ~ '^^/327rF^ [6], fiaton decay may produce too many 
(unthermalized) axions and upset standard nucleosynthesis. In the scenario of the fiatino 
LSP under consideration, it is also required that the LSP is not overproduced from fiaton 
decay. In general, a fiaton (p can decay into ordinary particles and their superpartners X, 
the axion a and the fiatino LSP Fi. After fiaton decay, the energy densities are given by 
Px = Bxp^ = (7r^/30)5fRj/T^H) Pa = BaP^ and pp^ = Bp^{mpjm^)p^ where p^ is the en- 
ergy density of the fiatons before they decay and the S's are essentially branching ratios of 
the fiaton decay. The precise relation between Ba and the fiaton branching ratios is given 
in [22]. As is well-known, nucleosynthesis (NS) puts a constraint on the extra amount of 
relativistic energy density, which is conveniently given in terms of the equivalently number 
of neutrino species, 6N^. At present [25], the bound is something like 6Ni, < 0.3 in the 
favored 'low deuterium' scenario. This can be translated to upper limits on Ba and Bp,^. 
Applying (p^/ pi,) < SNi, for the energy densities at the time of nucleosynthesis, one finds 
[19, 22] 

V / \ 1/3 

Ba<-^] SN^. (4) 
43 V43/4y 



3 



In the similar way, the energy density of the stable and massive LSP Fi has to be constrained; 
{PfiI Pi')ns < SNi,, which leads to 



7/4 my Tns ^ 



B,.<^^T^5N.. (5) 



To get a rough estimate of the reheat temperature, let us take the rate of the decay 
(/? — > X to be Fx ~ Fa/ Ba resulting in 

r.„ « - lOGev (^^) (a) 

where qrh ~ 100 is the effective number of particle degrees at Trh and we took the numerical 
factor Ba = 0.1. 

For Trh ~ 10 GeV, we get Bp^ < 10"*^ with Ba ^ 0.1. A more stringent limit on Bp^ 
comes from the flatino contribution to the present energy density. That is, in order to avoid 
the overclosure by the LSP, {pp^/pc)o < 1, one gets 

4^.^^ (7) 



Taking Trh ~ 10 GeV, one obtains Bp_^ ^ 10^^°. This is a very strong constraint. It appears 
impossible to get such a small number without relying on a severe fine-tuning of parameters 
in the model once flaton decay into flatino is allowed kinematically. Therefore, we have to 
forbid the decay of flatons into the flatino LSP imposing 

<2mp^. (8) 

Of course, the condition (7) can be invalidated in the case that the flatino is heavy enough 
to decay into ordinary particles, e.g., a neutralino and a light Higgs boson [22] which is 
opposite to our consideration. 

Forbidding the decay <^ — > Fi kinematically, we now have to fulfill the condition (4). 
Important thermalizable decay modes of a flaton to achieve Ba <C Bx ~ 1 include the 
decays into two top quarks, two light stops, two light Higgs bosons h^, two Z /W gauge 
bosons and a Z boson plus a Higgs boson. In this paper, we take the flatons to be as light as 
possible, so as to allow the lightest possible masses for the flatino LSP and for the NLSP. We 
therefore assume that only the last two decay modes are kinematically allowed. The decay 
mode into Higgs bosons comes from direct couplings of Higgs bosons and flatons given by 
Eq. (1) and has been considered in Ref. [22] in the context of an ordinary neutralino LSP. 
The other modes hsted above come from mixing between Higgs bosons and flatons induced 
by the same /i term interaction (1), which is the bosonic counterparts of mixing between 
ncutralinos and flatinos which has been worked out in Ref. [24]. (The qualitative features 
of the decay modes into top quarks or stops have also been considered in Ref. [19].) 

Under the condition of no CP violation, scalar flatons (denoted by F) can decay into 
two gauge bosons, e.g., F WW, and pseudoscalars (denoted by F') into a Z boson and a 
scalar Higgs, F' hPZ, as we describe below. The former mode gives a lower limit on the 
scalar flaton masses, 2mw < mp. Combining it with Eq. (8) we find 

2'mw < mp < 2m p^ < 2m^o (9) 
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where Xi is the NLSP neutrahno. (When mp < 2mw for instance, flatons can of course decay 
into three ordinary hght fermions mediating a sfermion. But the corresponding decay rates 
have a large phase space suppression which makes it hard to dominate over the decay mode 
F aa.) Concerning pseudoscalar flatons, there are more possibilities. When m^' > mp, 
the decay mode F' — > aF is open and has to be suppressed against the mode F' — > h'^Z. 
This requires 

m/jO + rriz < rrip/ < 2m^^ < 2m^o . (10) 

On the other hand, if mi?/ < mp, F' has no two body decay mode into an axion, and thus F' 
can be so light as to have the decay modes into light fermions; F' ff{Z) or ffh^. In this 
case, the corresponding decay rate has a large Yukawa-coupling or phase-space suppression 
factor leading to lower reheat temperature than in Eq. (6), and only the bound (9) is applied. 

Without resorting to a specific model, let us describe how the flaton couplings to gauge 
bosons arise. If we assume no CP violation, scalar flatons mix with the two CP-even Higgs 
bosons h^, and pseudoscalar flatons mix with the CP-odd Higgs boson A, through the 
/x-term potential as in Eq. (1). Quantifying small mixtures of flatons in Higgs bosons by 
^h°F, ^H°F s-nd Saf' in a self-explaining notation, it is straightforward to write down the 
flaton-Higgs couplings from the MSSM Lagrangian: 

L = -gmwevvF[FW+Wl^ + ^FZ^Z^] 

-— £^F'[cos(a - + sin(a - (3)H^%F'Z'' (11) 

Cw 

where Evvf = sin(ci; — (3)ehOp — cos{a — /3)eH0p, ol is the diagonalization angle of CP-even 
Higgses, 13 is defined by tan j3 — {H2) / (Hi) and cw — cos 9w with 9w being the weak mixing 
angle. 

Let us now consider the flatino-neutralino mixing [24] which arises at tree level due to 
the II term (1). For typical parameters, this mixing will give the main contribution to the 
flatino interactions with ordinary particles and their supcrpartners, dominating the one-loop 
processes [7] which are otherwise responsible for these interactions. The relevant interactions 
for our purpose are the decay of the NLSP to the flatino LSP, and the decays of sfermions 
to the flatino LSP. 

The term (1) leads to the mass matrix that mixes and {B, W3, Hi, H2) as follows; 

(0 iiSfidi iicp5i) , (12) 

where cp = cos j3, sp = sin f3, and the coefficients 5i ^ v/F^, depend on the specific form of the 
superpotential (1). Here v = 264 GeV is the Higgs vev. Let N be the diagonalization matrix 
of the MSSM neutralino mass matrix M^: Mj^ = N'^M^^N where Mjv is the diagonalized 
neutralino mass matrix for the eigenstate x°. Then, further diagonalizing (12) can be done 
by the small mixing elements; 

which describe the small content of the flatino Fj in the neutralino Note that we can 
have an enhancement in the mixing (13), s^qp. ^ lOv/Fa, when m^^p — mp. — 0(10) GeV. 
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Flaton-Higgs mixing parameters and flaton/flatino masses are dependent on specific forms 
of tlie terms (1) and (3). As we discussed, in order for tfie fiatino LSP scenario to be 
consistent witli cosmological considerations, tlie model sliould fulfill the mass relations (8), 
(9) and (10), and have sufficiently large mixing elements e in Eq. (11) to suppress the axion 
decay modes. However, note that mixing between the flaton sector and the MSSM sector 
is determined by the term (1) and flaton/flatino mass spectrum by the term (3) and its 
soft-breaking term. Therefore, arrangement for the mass relations (8,9,10) and large mixing 
elements e can be done by independent parameters in two separate sectors, which implies 
that such an arrangement can be easily achieved in generic flaton models. 

As the fiatino LSP has a mass in the 100 GeV region, its shght regeneration after the 
thermal inflation can provide a sizable contribution to the matter density at present. There 
are two important sources of the fiatino regeneration. The first is the neutralino decay [7] 
after its decoupling from the thermal bath. In the flaton models, the NLSP neutralino can 
decay into Z{h^)Fi or ffFi within 0.01 sec [24] and thus without affecting nucleosynthesis. 
Then, the ratio of the fiatino LSP density to the critical density is given by 

fi,. = ^Si,.. (14) 

For our case, mp^/m^o^ = 0{1) and thus we can have ~ fl^o ~ 0.1 — 1 which is 
valid in a wide range of the MSSM parameter space [26]. For this mechanism to work, 
the reheat temperature (6) after thermal inflation should be larger than the neutralino 
decoupling temperature ~ m^o/20. As can be seen from Eq. (6), this condition holds for 
low i%i ^ 10^° — 10^-*^ GeV which increases collider signals as well. 

A potentially more important source is thermal regeneration [27]. Since our reheat tem- 
perature is typically below 100 GeV, decay processes dominate over scattering processes for 
the thermal regeneration. To get a qualitative calculation of the fiatino population, let us 
take a sfermion-fermion-flatino interaction arising from the neutralino-flatino mixing, 

C = -geffpJfP, + h.c. (15) 

where Sffp^ oc v/Fg^ contains the factors from gauge quantum numbers and neutralino- 
flatino mixing and /// denotes a left-handed or right-handed sf ermion/f ermion field. The 
thermally regenerated fiatino population is given by VLp^h? = mpYp^/'i.bb eV with the factor 
[27, 28] : 



p^2x 10"^ — 



-RH 

-,2^2 



Fi ~ - - -7^F{xj) (16) 



where F^ = g e jj^p^mj/Sn is the decay rate and F{x^) is the Boltzmann suppression factor 
as a function of x^ — m^/T^Yi given by 

Fixf) = \ I ^ I - - tan"^ , ^ 1 x^ + x? (x^ - 2XfJx^ - x'^,) 
^ x^^ Jxf - 1 [\2 - x)) f) 

We have defined the function F{xj) to remove the dependence on the axion scale Fa for the 
remaining factor; MpV ^/T^-^ oc ^ f/^^p where F^ is the flaton decay rate. To a very good 
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approximation for the range of a; ^ 5, the function F(x) can be expressed in terms of an 
analytic function; F{x) ~ 36e~°'^^^ giving 

a,_ft^«5xl0'(^)(ip)'4e-»-/. (17) 

From the above equation, one finds that the requirement Vtp^h^ ^ 1 gives Xj^l7 with m^^ = 
and ffjjp^ = 10~^. This consideration gives some meaningful bounds on sfermions masses, 
TTif > 17Trh given the reheat temperature (6), or vice versa. The thermally regenerated 
population of the fiatino LSP can also be the cold dark matter which, however, needs a 
fine arrangement for the sfermion mass and reheat temperature as VL is a very sensitive 
function of xj. Note that the thermal regeneration becomes easily negligible as Fa becomes 
large. For instance, for Fa = 10^^ GeV, one has Trh ~ 2 GeV (6) and thus > 50 for 
ruf > 100 GeV. 



III. A MODEL 

To illustrate our discussion more explicitly, we consider the minimal case of two fields 
and n = 4, which has been analyzed in Ref. [22]. This model contains two scalar flatons, 
one pseudoscalar flaton, and two fiatinos. The flaton superpotential is 

W,^ = ^/'Q (18) 

with U{l)pQ charges of P and Q being 1 and —3, respectively. Analyzing the scalar potential 
including soft supersymmetry breaking terms 

Vsoft = ^P'Q + h.c., (19) 

one finds scalar flatons F2^i with masses-squared 

= ^/Vo (3(12 -0 + ^'(12 + ± |12 - C|Va;4 + 42a;2 + 9) (20) 
and a pseudoscalar flaton F' with mass-squared 

ml,^ffilax' + 9) (21) 
where the parameters are defined by 



In the above, the flatons P, Q are expanded as 



P=-L((P) + P'-3.fF') 

Q = 71 (^^) + ^'-^^^') (22) 
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where = (P)^ + 9((5)^. The diagonahzation matrix of the scalar flatons P' and Q' is 

(23) 

where the mixing angle </? is determined by 



P' \ _ / cos (/? — sin (/? \ / F2 
Q'y vsin(/? cos(/? ) \Fi 



+ 3 

cos2(^ = S5in(12 — ^)- 



Va:;"^ + 42^2 + 9 



6x 

si„2^ = .,n(12-0-^__. (24) 



Prom the superpotential (18), one also finds the fiatino masses 



SffioiVx^ + 1 ± 1] (25) 



The light flatino Fi is supposed to be the LSP. The rotation matrix from the flavor states 
P, Q to the mass eigenstates F2, Fi is given by 



P\ _ f cos ip — siTnp\ / Fi 
Q J V sin cos <^ y V -^2 



(26) 



with the mixing angle (f satisfying 



1 X 
cos2(^= — ^==, sm2(f^ — TTTt- ^^^^ 

V^c^ + 1 V 2; + 1 

As all the masses contain the overall factor |//io|, it is now straightforward to flnd the 
region of the parameters x and ^ satisfying the desired mass relations. It turns out that 
mFi2-i 'iTT'F' < '^TTip^ requires 

(A) e<12; x>2.6, Ub<^<^c 

{B) i>12- X > 3.5, 12 < e < Co (28) 
VFT42x2T9 - - 3 



with = 12 

^6 = 12 

^0 = 36 



Va;^ + 42x2 + 9 + ^2-3 

+ 42^2 + 9 + 12Va;^ + 1 - 5^2 
Va:;^ + 42x2 + 9 + ^2 - 3 
a;2 + 2 - 2^/^^ 
^2 + 9 



We always have mp^ < mp', mp^ in the regions (A,B) and also mp^ < mpi in the region (B). 
The requirement (10) has to be fulfilled in both regions. 

The important interactions of fiatons (fiatinos) and MSSM fields comes from the mixings 
between flatons (flatinos) and Higgs bosons (neutralinos) due to the /x-term superpotential 
of our choice 

PO 

h-^H,H2 (29) 

IVl p 

with U{l)pQ charge of H1H2 being +2. Note that fi = h{PQ)/2Mp = hfj^Q. First of aU, 
Eq. (29) determines the flatino-Higgsino mixing elements 6i deflned in Eq. (12) as follows: 



()i = -7; {c^ + xs^}, 02^— {-s^ + xc^), (30) 

-Ta, X -Ta 
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with = cos,^ and = sin,^. Including the soft term of (29), we get the Higgs-flaton scalar 
potential, 



h 



PQ 



+ |i/2r + 



h 



PQ 



M2 



+ c.c. 



We note that the CP-odd Higgs boson has mass-squared 



(31) 



sin 2/3 \ fj, h 
It is then straightforward to find the following fiaton-Higgs mixing: 



{— sin(Q; — /5)[2s^ — 2xc^ 



A, 



-Fcos(q; - /3)[-(4x + 6s<^ - Qxc^) + — {s^- xc^)]} 



1^ 



^H^Fi — 



V 



„ „-{cos(q; + (5) [2s^ - 2xcJ 
X mjjo — rrip^ ^ ^ 

-Fsin(Q; (3)[j{4:X^s^ -\- 6s<^ - 6xc^) + -^{s^- xc^)\} 



X mlo - mj,^ 
J 



{sin(Q; — P)[2c^ + 2xs^ 



A, 



- cos{a - (3)[-{4:X + 6c^ + 6xs^) + — (c^ + xs^)]} 



1^ 



X m 



-{ - cos(q; + /3)[2c^ + 2xs^ 
A, 



sm{a + /3)[-{4:X + 6c^ + 6xs^) H (c<^ + xs^)]} 

lb jJi 



f x' + 3 jj? 



X m\ — rrip, jj, 



(Ah / , 

2 {—-67:} 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



where v = yJ{Hiy + {H2y. The flaton decay rates for the processes F — > WW and F' 
hPZ of interest are given by 



V{Fi WW) = a2\evvF, 



1 - 4- 



m 



w 



m 



(38) 



Fi 



0.2 



cos(q; — (5Y\eAF'\^'mF' [ 1 + 



2 12^ 

^ftO ~ 2^^ 

ml, 
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where SyvFi = sin(Q; — (3)ehop. — cos{a — (5)eH0Fi- The above rates have to be compared with 
the rates for the decay processes F ^ aa and F' aF^ [restricting ourselves to the region 
(B) in Eq. (28)] which are given by 



r(Fi aa) 



1 m-j,^ {-xs^ + 9c^) 
32^"Ff (a;2 + 9) 



\2 



mi (xa„ + 9s„ 



\2 



r(F' ^ aFi) = 



1 m|„ / (3c^ - 3a;s^)^ 



167r F2 \^ ml, J {x^ + 9) 



Let us define R{Fi) = T{Fi aa)/r{Fi WW) and R{F') = r{F' aFi)/V{F' h^Z) 
which behave hke 



R{F,)o,^'^, R{F') oc ^ ~f ^'^ (40) 



apart from the other numerical factors. Eq. (40) shows that small R can be easily obtained 
with small fiaton masses and large /j,. In the case of small /i, one would need to have small 
rriA as well as mi?/ < v. 

Having calculated all of the relevant quantities, we could in principle identify the region 
of parameter space which gives a viable cosmology as well as viable particle physics. Here, 
we present instead two sets of parameters, with (i) a bino-like and (ii) Higgsino-like NLSP, 
which satisfy all the requirements. The sets give an NLSP mass around the lowest possible 
value (m/jo +mz)/2. The flatino LSP is only a little fighter, while the other flatino as well as 
the flatons have masses around 200 GeV, features which appear to be typical for the region 
of parameter space corresponding to a light NLSP. 

In making the choice of parameters, we focussed on the region 2m]y < rriFi, < 2m2, 
which can be arranged for ,^ ~ 12 [see Eq. (20)]. Both parameter sets have 

tan /? = 3 , rriho = llOGeV , M2 = 2Mi , 

and the other parameters are listed in in Table I along witfi tfie important output quantities. 
In both cases, we have R{Fi, F') < 0(0.01) and the reheat temperature Trh = 16 GeV using 
the decay rates (38). The neutralino NLSP decouples after refieating, and can decay into 
the flatino LSP to provide the cold dark matter of the universe. To check that potential 
overproduction of the flatino LSP from the decay of thermal sfermions can be avoided, let 
us take as an example the right-handed stau. From e^op-^ in Table I, the coupling in Eq. (17) 
is given by £f^T-B.h ~ ^'"^ ^ wfiicfi gives Qp_^h'^ ^ 1 for mf^ ^ 300 GeV for the case (i) 
and a similar figure is obtained for the case (ii). 



IV. CONCLUSION 

In an interesting class of models, a non-renormalizable term of the superpotential is 

responsible for the spontaneous breaking of PQ symmetry, while another such term generates 
the fi term of the MSSM. The fiaton fields which break PQ symmetry are accompanied by 
fiatinos, and the lightest fiatino can be the LSP. Through the term of the superpotential, 
the NLSP neutralino decay might then be visible at coUiders. 



10 



In this paper, we have examined the cosmology of this kind of model, following our earlier 
work [22] which explored the same flaton models on the assumption that all of the fiatinos 
were unstable. We make the reasonable assumption that the fiaton fields are nonzero in 
the early Universe, giving thermal inflation which eliminates pre-existing relics and leads to 
a rather well-defined cosmology. The decay of flatons into relativistic axions will interfere 
with nucleosynthesis unless the branching ratio is small, while the decay of flatons into 
flatinos should be forbidden altogether so as to avoid overproduction of the flatino LSP. We 
have argued that these combined requirements lead to the bounds m^o > mp^ > mw or 

{mz + m/io)/2 on the masses of the neutralino NLSP Xi ^iid the flatino LSP Fi. In order 
to have the hghtest possible NLSP, one can focus on the case of fight scalar (pseudoscalar) 
flatons, which can decay only into two W bosons (a light Higgs and a Z boson). These 
decays, as well as the neutralino NLSP decay into the flatino LSP, come from the mixing 
between flatons (flatinos) and Higgses (Higgsinos) due to the n term of the superpotential. 

We focus on the case of a low axion scale. Fa ~ lO^*^ GeV, which increases the feasibility 
of observing NLSP decays in future colliders. While the axion can hardly be the dark matter 
with such a low scale, the flatino LSP becomes a good dark matter candidate, because the 
reheat temperature can be high enough to generate the NLSP in thermal equilibrium, but 
low enough to suppress fiatino production from the decay of thermally produced sfermions. 

To make some definite statements, we studied an explicit model with F^ ~ 10^*^ GeV. 
We verified that flaton decays into WW or h^Z can be large enough compared to their 
decays into axions in a wide range of parameter space, giving rise to the reheat temperature 
Trh ~ 0(10) GeV. In this model, the NLSP can be as fight as ~ 120 GeV with m,,o = 110 
GeV, while the other fiatino and the fiatons have masses ^ 200 GeV. We further notice that 
the parameter space admits the interesting mass region, m^o — mp^ = 0(10) GeV, where 
an enhancement of the fiaton-neutralino mixing, s^op^ ~ 10~^, results and thus the NLSP 
decays becomes faster. 
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Biuo-likc XLSP 


niggsiuo-likc XLSP 


Inputs 


Ml 


120 GcV 


220 GeV 




287 GcV 


130 GeV 


TUA 


300 GeV 


150 GcV 




305 GeV 


162 GcV 


X 


4 


4 


f/h 


1/24 


1/11 




-1 


-2.2 




13/24 


13/11 


Outputs 




6.6 X 10^° GcV 


4.4 X 10^° GeV 


niF' 


216 GeV 


223 GeV 




175 GeV 


181 GcV 




162 GeV 


168 GcV 




184 GeV 


190 GeV 




112 GeV 


116 GeV 


m^o 

Xi 


121 GeV 


123 GeV 




-6.4 X 10-^ 


37 X 10-8 




(15,-1.7,15.6,-4.8)?;/Fa 


(10,-6.1,-44,-41)t;/Fa 




16 GeV 


16 GeV 



TABLE L Two representative parameter sets taking tan/? = 3, m^o = 110 GeV, = 2Mi and 
^ = 13. Note that x° = {B, W3, -ff? , H^) and v = 264 GeV 
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